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Abstract 

Research of delayed neural networks with variable self-inhibitions, inter- 
connection weights, and inputs is an important issue. In this paper, we discuss 
a large class of delayed dynamical systems with almost periodic self-inhibitions, 
inter-connection weights, and inputs. This model is universal and includes de- 
layed systems with time- varying delays, distributed delays as well as combi- 
nation of both. We prove that under some mild conditions, the system has a 
unique almost periodic solution, which is globally exponentially stable. We pro- 
pose a new approach, which is independent of existing theory concerning with 
existence of almost periodic solution for dynamical systems. 

Key words: Delayed Dynamical Systems, Almost Periodic Solution, Global 
exponential convergence. 

1 Introduction 

Recurrently connected neural networks, sometimes called Grossberg-Hopfield neural 
networks, have been extensively studied in past years and found many applications in 
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different areas. However, many applications heavily depend on the dynamic behaviors 
of the networks. Therefore, analysis of these dynamic behaviors is a necessary step 
toward practical design of these neural networks. 

Recurrently connected neural network is described by the following differential 
equations: 

—JT- = -diUi(t) + J2 a ij9j(uj(t)) + h i = l,---,n (1) 

j 

where gj(x) are activation functions, d^ are constants and Ii are constant inputs. 
In practice, however, the interconnections are generally asynchronous. Therefore, one 
often needs to investigate the following delayed dynamical systems: 
duAt) n n 

— l — = -d iUi {t) + Y, a ij9j( u j(t)) + S kjfiMt ~ T ij(t))) + k i = 1, • • • ,n (2) 

0,1 3=1 3=1 

where activation functions gj and fj satisfy certain defining conditions, and a^, b^, Ii 
are constants. 

However, the interconnection weights a^-, bij, self-inhibition di and inputs Ii should 
be variable with time. Therefore, we need to discuss following dynamical systems. 
du- 11 n 

— - = -di(t) Ui (t) + Y aaWgjMt)) + E h.jlDfjioAl - Tij(t))) + Ii(t) (3) 

0,1 3=1 3=1 

or delayed systems with distributed delays and time-varying delays 
du(t) n 

— ^ = -di(t)ui(t) + 

ai j=1 

n PCX) 

+ J2 h iA l ) / kMfjM* - T ij(t) - s))ds + Ii(t) % = 1, • • • , n (4) 

3=1 ° 

To unify these delayed systems, in [2], we proposed the following general dynamical 
systems 

dui n 

-± = -di(t) Ui (t) +Y a ij( t )9j( u j( t )) 

+ E / fiM* ~ T i 3 \t) - s))dKij(t, s) + Ii(t) % = 1, • • • , n (5) 

.7 = 1 ° 
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where > and s — > dKij(t, s) is a Lebesgue-Stieljies measure for each t E R. 
The initial conditions are 

Ui(s) = <fii(s) for se (-oo,0] (6) 

where (pi G C((— oo, 0]), i — 1, • • • , n. 

There are several papers discussing periodic dynamical systems and their periodic 
solutions and its stability, see [1, 2, 6, 4, 7, 8, 9, 10] and others. For example, in [9, 10], 
under assumptions that di(t) > d t > 0, ciij(t),bij(t), Ii(t) : R + — > R are continuously 
periodic functions with period uj > 0, i.e., di(t + u) = di(t), = a^(t + a;), 

bij(t) = bij(t + u), Ii(t) = Ii(t + u), Tij(t + uj) = Tij(t) for alH > and i, j = 1, 2, . . . , n, 
the model (2) was investigated. Existence of periodic solution and its stability were 
investigated. And in [9], comparisons of various stability criteria were given too. In 
[2, 3], a general approach to discuss existence of periodic solution and its stability 
analysis for the model (3) was proposed. In [11, 12], authors presented some results 
on almost periodic trajectory and its local attractivity of shunting inhibitory cellular 
neural networks (CNNs) with delays. In [13], authors proved existence and attractivty 
of almost periodic solution for CNNs with distributed delays and variable coefficients. 

In this paper, we will investigate the existence of almost periodic solution of (5) 
and its global stability. All the conclusions and methods used apply to discussion of 
periodic solution when d iy bij etc. are periodic functions or equilibrium point when 
di, dij, b^ etc. are constants. 



2 Some preliminaries 

Definition 1 Class H{Gi, • ■ ■ , G n } of functions: Let G = diag[Gi, ■ ■ ■ , G n ], where 
Gi > 0, i = 1, - ■ ■ ,n. g(x) = (gi(x), ■ ■ ■ , g n (x)) T is said to belong to H{d, • • ■ , G n }, 
if the non-decreasing functions i = l,---,n satisfy !M£±B^zMgH < Q. Also, 
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Hi{Fi, • • • , F n } of functions: Let F = diag[Fi, ■ • • , F n ], where Fj > 0, i — 1, • • • , n. 
f(x) = (/i(x), ■ • • , f n (x)) T is said to belong to • • • , F n }, if the functions fi(x), 



Definition 2 Let x(t) : i? — > R n be continuous, is said to be almost periodic 
(Bohr) on R if for any e > 0, it is possible to find a real number I = /(e) > 0, for any 
interval with length /(e), there exists a number uj = u(e) in this interval such that 
\\x(t + uj) -x(t)\\ < e, for Vt G R. 

Throughout the paper, we use following norm 



In this section, we investigate the dynamical system (5) with initial condition (6). In 
the sequel, we make the following assumption. 

Assumption A 

1. Activation functions g(-) G H{G ± , G 2 , ■ ■ ■ , G n }, /(•) G H 1 {F 1 , F 2 , ■ ■ ■ , F n } 

2. di(t), a,ij(t), Tij(t) > and Ii(t) are continuous, and di(t) > di > 0, for i,j = 



3. For any s G R, Kij(t,s) : t — > Kij(t,s) is continuous, and for any t G R, 
dKij(t, s) : s — > Kij(t, s) is a Lebesgue-Stieljies measure, for all i,j — 1, 2, • • • ,n 

4. For any e > 0, there exists / = /(e) > 0, such that every interval [a, a+l] contains 
at least one number u> for which 




U (t)\\{(,oo} = .™ aX J& VWI 



3 Main Results 



n 



diit + u) -di(t)\ < e 
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\aij(t + uj) - a,ij(t)\ < e 
\Ii{t + u) - Ii(t)\ < e 

\Tij(t +Uj) - Tij(t)\ < e 

/ \dKij(t + co,s) — dKij(t, s) | < e 

J 

hold for all i, j = 1.2. • • • , n and t E R. 

5. \dKij(t,s)\ < \dKij(s)\, and for some e > and J °° e es \dKij(s) \ < oo . 

It can be seen that under item 4 in the assumption A, di(t), ciij(t), Ii(t) and Tij(t) are 
almost periodic functions. Therefore, they are all bounded. We also denote |a*J = 
su P{teR} K0)| < °o, \ b *ij\ = sup {tefl} / °° \dK i:i (t,s)\ < oo, \I*\ = sup {teR} \Ii(t)\ < oo, 
t*j = sup {teR} Tij(t) <oo, i,j = !,-■■ ,n. 

Lemma 1 Suppose that all items in Assumption A are satisfied. If there exist 
£i > ; i — 1, • • • , n, such that 

n n r-OG 

-dim + E K(t)\G^ + £ F& / \dK ZJ (t, s)\ <- v <0 (7) 
/io/d /or all t > 0. In particular, if 

n n poo 

-d& + E K-|g&- + E / < ( 8 ) 

TTien any solution u(t) of the system (5) is bounded. 

Proof: Define M(t) = max ||m(s)||{£,oo}- It is obvious that ||{^,oo} < M(t), and 
M{t) is non-decreasing. We will prove that M(t) < max{M(0), |J}, where 

7 = max{j/*| + £ [|a^j|^(0)| + l^-j |/,(0)|1 } 
For any fixed to, there are two possible cases. 
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Case 1. 



\u 



(*o)||{$,oo} < M(t ) = max||w(s)|| {?i0o} 



s<to 



In this case, in a small neighborhood (t ,t + 5) of t , \\u(t) ||{^ )00 } < M(t ), 
M(t) = M(t ) 



Case 2. 



IM*o)||tf,oo} = M(to) = max ||u(s)|| tti0o} 

S^Zq 



In this case, let i to be such an index that 

CoKo^o)! = ||M(*o)|ltt,oo} 

Then noticing 

\9j{*)\ < Gi\*\ + \9i(P)\ \fi(*)\<Fi\*\ + \fi(0)\ for j = l,---, nseR 
we have 

r n 

= sign(u itQ (t Q )) - d ito (t )u itQ (t ) + £ ^3^) 9 3^3%)) 

1 3=1 

n roo 

+ £ / /i(«j(*o - 7itoi(*o) - s))dK itQj (t , s) + 4 o (t ) 
i=i J 

n 

+ / l^(io-T^(to)-s)|e-V^(^,s)| + l4 (^)l 



J'=l 



+ El a *w(*)IM°)l + / |d*w(*.«)ll/i(0)l 



J'=l 



< 



-4>)6 t0 + £l<JG,< 



+ E%7 l^ i(*o>s)l ll«(*o)|ltt,oo} + / 
j"=i J 

< -n\\u{t Q )\\ {too} + i = -nM{t ) + i 



Thus, if M(to) > ||/|, then M(t) is decreasing in a small neighborhood (t ,t + S±) of 
t - On the other hand, if M(t ) < -\I\, then in a small neighborhood (t ,t + 5 2 ) of 
t , ||^(^)||{?,oo} < Therefore, M(t) < max{M(t ), ^|/|} in a small neighborhood 
(t , t + 5) of t , where 5 = min{6i, <5 2 ). 

In either case, we have M(t) < max{M(t ), ||/|} in a small neighborhood (t , t +5) 
of t . 

Starting from t = 0, process this argument at every t > 0, we conclude that 
if M(0) > ||/|, then M(t) = M(0) for all t > 0. Instead, if M(0) < ^|J|, then 
M(t) < for all t > 0. Therefore, M(t) < max{M(0), -|/|} for all t > 0, which 
proves that u(t) is bounded. Lemma 1 is proved. 

Lemma 2 Suppose that all items in Assumption A are satisfied. If there exist 
£i > ; i — 1, 2, • • • , n, and (3 > 0, such that jo all t > 0, there hold 

n n roc 

-diim + E \a tJ (t)\G^ + E / e^|^(t, s)| < -r, < (13) 

T/ien /or any e > 0, t/iere exist T > and I = 1(e) > 0, such that every interval 
[a, a + I] contains at least one number uj for which the solution u(t) of system (5) 
satisfies 

\\u(t + u) - u(t)\\{( i0O } < e forallt>T (14) 

Proof: Define 

n 

€i(u,t) = -[di(t + di(t)]ui(t + uj) + J2[aij(t + u) - aij(t)]gj(uj(t + u)) 

3=1 

n roo 

+ E / - ^A 1 + s))- fj(uj(t - Tij(t) +lu- s))]dK ij (t + u, s) 

j=l J0 

n roo 

+ E / fjM* - T ij(t) + u - s))d[Kij(t + w,s)- Kij(t, s)] + [hit + uj)- h(t)} 
j=l J0 

By Lemma 1, u(t) is bounded. Thus, the right side of (5) is also bounded, which 
implies that u(t) is uniformly continuous. Therefore, by the assumption A, for any 
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e > 0, there exists I = 1(e) > such that every interval [a, a + 1], a e R, contains an 
uj for which 

(15) 



\ei(u, t)\ < —rje for all t E R , i — 1, 2, • • • , n 



Denote Xi(t) = Ui(t + u) — Ui(t). We have 
dxi(t) 



dt 



-di(t)xi(t) + + 9j(uj(t))} 

3=1 



n POO 

+ E / [foMt + u - T ij(t) - s )) - fj( u j(t - T ij(t) - s))]dKij(t, s) 

+6i(u,t) (16) 

Let i t be such an index that 

tu\xH(t)\ = Mt)h,oo } (17) 

Differentiate e^ s |x it (s)|, we have 

Pe^\x it (t)\ + e^sign(x it (t)){ - d k (t)x it (t) 



s=t 



+ E 9j(uj(t + w)) - gj(uj(t)) 

3=1 1 



n rc 



fj(Uj(t + U- T itj (t) -S))- fj{Uj(t ~ T itj (t) - s)) 



dK kj (t,s) + e it (u,t) 



< _ _ p]\x it mr t % + E kdWI^M*)^ 

1 3=1 



3=1 



+ E% If - r w(*) - a)ie7 1 c-^W + '^^J|(«ir w (t J a) | } + ^/ee* (18) 



Similar to the proof of Lemma 1, let 

*(*) = m fXp1z(s)ll{£,oo}} 



(19) 



If there is such a point to > that ^(io) = e/3 *° 11^(^0) ||{^,oo}- Then we have 

'V"K,oi; 

t=t 



(20) 



In addition, if ^(t ) > ee' 3 * , then \&(t) is decreasing in a small neighborhood (t , t + 8) 
of t . On the other hand, if ^(t ) < ee' 3 * , then in a small neighborhood (t ,t + 5) of 
to, e /3 *||a;(t)||{^ i00 } < ee' 3 * , and < max{^(t ),ee /3t °}. By the same reasonings used 
in the proof of Lemma 1, for all t > to, we have e**||x(i) ||{£,oo} < max{^(to), ee' 3 *}. 
Therefore, there exists ti > 0, for all t > ti, \\x(t)\\^ tOQ y < e. 

Instead, if for all t > 0, we have V(t) > e pt \\x(t) ||^ j00 }, then = ^(0) is a 
constant, and e /3 *||x(t)||{^ 00 } < ty(t) = ^(0). Hence, there exists t 2 > 0, for all t > t 2 , 
such that || ^(t) || {£,00} < e - 

In summary, there must exist T > such that ||^(t)||{^ j0 o} < e holds for all t > T. 
Lemma 2 is proved. 

Theorem 1 Suppose that all items in Assumption A are satisfied. If there exist 
£i > ; i — 1, 2, • • • , n, and f3 > suc/i that 

- [diW-^i + EM*)l^ + E%-^ / e^|rf^(t,s)|<0 (21) 
/ioW /or all t > 0. In particular, 

n n />oo 

- [* - + E + E / ^|^.( S )| < (22) 

T/ien the dynamical system (5) has a unique almost periodic solution v(t) = [v i(t), v 2 (t), 
■ ■ ■ ,v n {t)} T , and for any solution u(t) — [ui(t), u 2 (t), u n {t)] T of (5), there holds 



\u{t)-v{t)\\ = 0{e-*) (23) 



Proof: As did before, define 



Ci,fc(f) = ~[di(t + t k ) - di(t)]ui(t + t k ) + Et^i^ + *fc) - a>ij(t)]9jMt + t k )) 

j'=i 

+ E / [fj( u i(t ~ T ij(t + **) + *fc - s )) - - r u(0 + *fc - + h, s) 

, r'" 

n /-oo 

+ E / + ** - - s )M^(* + h, s) - K^t, s)\ + [h{t + * fc ) - h{t)\ 

(24) 



where tk is any sequence of real numbers. From the Assumption A and the bounded- 
ness of u(t), we can select a sequence {tk} — > oo such that 



h,k(t)\ < - for all i,t 



(25) 



Because {u(t + tk)} k x L 1 are uniformly bounded and equiuniformly continuous. By 
Arzala-Ascoli Lemma and diagonal selection principle, we can select a subsequence 
t kj of tk, such that u(t + £ fe ) (for convenience, we still denote by u(t + tk)) uniformly 
converges to a continuous function v(t) = [vi(t) , v 2 (t) , • • • , v n (t)] T on any compact set 
of R. 

Now, we prove v(t) is a solution of system (5). In fact, by Lebesgue dominant 
convergence theorem, for any t > and St G R, we have 

Vi(t + St) - Vi(t) = lim Ui(t + 5t + t k ) - Ui(t + t k )) 

/■t+5t r n 

= lim/ \ - di(a + t k )ui(a + t k ) + y2 a ij( a + t k)9j(u j (a + t k )) 

k^ooJt L j . =1 

TO /-OO "I 

+ / fj( u j( a + tk- nj(o- + ifc) - s))dKij(a + t k ,s) + Ii(a + t k ) >da 

3=1 ° ' 
rt+St r n 

= / - di(a)ui((T) + J] M^ffiM* 7 )) 
* 1 j"=i 

+ S / ~ r ii( a ) ~ s))dKij{a, s) + /;(<r) W<7 

i=i J 

+ lim / e iik (s)da 

/■t+<St r n 

= / - di(a)vi(a) + a i3i a )9j{vj{o-)) 
1 1 j=i 

n POO ~\ 

+ J2 fi( v A a - T iA°) - s))dK i:j (a, s) + h(a) \da 
3=1 ' 

which implies 

-di(t)vi(t) + J2 a i3(t)93Mt)) + r fjMt-T tJ (t)-s))dK l3 (t, s)+Ii(t) (27) 

.7=1 ° 



(26) 



dv., 
It 
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Therefore, v(t) is a solution of system (5). 

Secondly, we prove that v(t) is an almost periodic function. By Lemma 2, for any 
e > 0, there exist T > and I = /(e) > 0, such that every interval [a, a + I] contains 
at least one number oj for which 

\ Ui (t + u) - Ui(t) | < e for all t > T (28) 

Then we can find a sufficient large K e N such that for any k > K 

\ui(t + t k + u) - Ui(t + t k )\ < e forallt>0 (29) 

holds. Let k — > oo, we have 

\vi{t + u) -Vi(t)\ < e forallt>0 (30) 

In other words, v(t) is an almost periodic function. 

Finally, we prove every solution u(t) of (5) converges to v(t) exponentially with 
rate (3. 

Denote y(t) — u(t) — v(t), we have 
dv(t) n 

Let z t be such an index that 

\yiM\-^\\y(t)h,oo } (32) 

and differentiate e^ s \yi t (s)\, we have 

^{ e/3S |^( s )l}| s=i = ^IWtWI + e*«Wlfe(*)){ - d*(*)lfe(*) 

n 

+ £ Oi« (*) ( u j (t) ) - & fa (f ))] 

3=1 

11 



n s 

+ Y,ifj( u j( t - t <«j(*) - s )) - - T itj(t) - s))]dK itj {t, s) \ 

3=1 J 

k j=i 

+ / ly^-r^-W-^I^V^^W^^I^t.a)!} (33) 

j"=i ' 

Define A(t) = max|e /8s ||y(s)||{^ j00 }|. If for some t , 

A(*b)=e^||y(to)|| {€ ,oo} (34) 

Then, by (33), we obtain 

j t {e^\y U) (t)\}t=t <-V^(to)<0 (35) 
and A(t) — A(t ) in a small neighborhood (t , t + 8) of t . On the other hand, if 

e /3io |b(to)||«,oo}< A(t ) (36) 
Then in a small neighborhood of t , 

e*||y(t)|| {€>oo} <A(to) (37) 
Therefore, A(t) = A(0) for all * > 0, and 

Ib^ll^ < A(0)e"^ (38) 

Namely, 

IK*)-«(*)ll«,oo}< A(0)e^ (39) 

Theorem 1 is proved. 

Remark 1 From the properties of continuous and almost periodic functions, i.e., 
Bochner's Theorem, see (pp. 4-7, [5]) it can be seen that if di(t), %(£), hj(t), Tij(t) 
and are all continuous and almost periodic, then the item 4 in Assumption (A) 
are satisfied. 
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If dKij(t, 0) = bij(t), and dKij(t, s) = 0, for s ^ 0, then (5) reduces to 
du- n 11 

— - = -di(t)ui(t) + E a^gjMt)) + E bi^fjMt - Tij (t))) + h{t) (40) 

al 3=1 3=1 

Therefore, we have following 

Corollary 1 Suppose g(x) = (g^x),---, g n (x)) T G H{G U • • • , G n } and f(x) = 
(fi(x),---J n (x)) T G H{Fi, • • • , F n }. di(t), otjit), bij(t), Tij(t) > and I t (t), i,j = 
l,---,n, are continuous almost periodic functions. If there are positive constants 
£i, • • • , £ n and (3 such that 

n n 

U-di(t) + (3] + E^K-WIGj + E |^(t)|F,e^ < 0, i = 1, • • • ,n, (41) 

3=1 3=1 

hold for all t > 0. In particular, if 

n n 

+ /3]+E &KI G J + E l&yl V« < 0, i = 1, • • • ,n, (42) 

J=l 3=1 

Then the dynamical system (40) has a unique almost periodic solution v(t) and, for 
any solution u(t) of (40), there holds 

||«(t)-t;(t)|| = 0( e -*) (43) 

Instead, if dK^t, s) = bij(t)kij(s)ds, then we have the following distributed delayed 
dynamical systems: 

du(t) n 

— ^ = -di(t)ui(t) + 
ttC j=i 

+ E / M s )/i( u i(* - r *i0) - s )) rfs + J i (*). i = 1, • • • , n (44) 
j"=i 

Corollary 2 Suppose £/(:r) = (<7i(x), ■ • • , g n (x)) T G #{Ca, • • • , G n } and f(x) = 
(fi(x),---,f n (x)) T G H{Fi, • • • , F n }, di(t), aij (t), bij(t), hj(t), Tij (t) > and L t (t), 
i,j = 1, • • • ,n, are continuous almost periodic functions and there exists (3 > such 
that 

J e ps \k i:j (s)\ds < oo (45) 
13 



If there are positive constants £1, • • • , £ n such that 

n n POO 

a-di + 13] + E&K'WI^ + E M*)!*^ / e*|M 5 )|ds < °> ( 46 ) 
i=i i=i 

hold for alH > 0. In particular, if 

n n />oo 

+ + E&KI G J + E I ^iV* / e^|fey( a )|da < 0, (47) 

j=l 3=1 ° 

Then the dynamical system (44) has a unique almost periodic solution v(t) and, for 
any solution u(t) of (44), we have 

|Kt)-^)|| = 0(e^) (48) 

Remark 2 In [11, 12, 13], authors discussed almost periodic solution for CNNs 
with delays (constant delays in [11], = r(t) in [12], and distributed delay with 
Tjj(t) = in [13]) and its local stability. Our model unifies models with time- varying 
delays Ty(i) and distributed delays and even both, and includes all the systems in these 
three papers as special cases. The stability is global. Moreover, in these three papers, 
all the reasonings heavily depend on the Theorem 1.93 in [14]. Thus, many unrea- 
sonable restrictions are imposed on the activation functions and the interconnection 
weights (for example, interconnection weights must satisfy [%•(£) | > c and | &«•,•(£)[ > c 
for some positive constant and many others). In this paper, we propose a totally dif- 
ferent approach to prove existence of almost periodic solution and its global stability 
without citing any complicated theory. As we show, all the redundant restrictions 
needed in these papers are removed. 

References 

[1] J. Cao, "New results Concerning exponential stability and Periodic solutions of 
delayed cellular neural networks" , Phys. Lett. A, 307, 136-147, 2003 



14 



[2] T. Chen, W. Lu and G. Ron Chen, "Dynamical behaviors of a large class of 
general delayed neural networks" , Neural Computation (2005) in press 

[3] W. Lu and T. Chen, "On periodic dynamical systems", Chin. Ann. Math., 25B:4, 
455-462, 2004 

[4] K. Gopalsamy and Sariyasa. "Time delays and stimulus-dependent pattern forma- 
tion in periodic envirorments in isolated neurons", IEEE Trans. Neural Networks, 
13:2, 551-563, 2002 

[5] B.M. Levitan, and V.V. Zhikov, "Alomost periodic functions and differential 
equations", Cambriadge University Press, New York, 1982. 

[6] S. Mohamad, and K. Gopalsamy, "Neuronal dynamics in time varying environ- 
ments: continuous and discrete time models", Dis. Cont. Dyn. Syst, 6, 841-860, 
2000, 

[7] L. Xiang, J. Zhou, Z. Liu, "On the asymptotic behavior of Hopfield neural network 
with periodic inputs", Appl Math and Mech, 23:12, 1367-1373, 2002 

[8] Zhang, Y. "Absolute periodicity and absolute stability of delayed neural net- 
works", IEEE Trans. Cirsmts and Syst. I, 49, 256-261, 2002 

[9] Y. Zheng and T. Chen, "Global Exponential Stability of Delayed Periodic Dy- 
namical Systems", Phys. Lett. A, 322:5-6 344-355, 2004 

[10] J. Zhou, Z. Liu and G. Ron Chen, "Dynamicas of delayed periodic neural net- 
works", Neural Networks, 17:1, 87-101, 2004 

[11] X. Huang and J. Cao, "Almost periodic solution of shunting inhibitory cellular 
neural networks with time- varying delays" ,Phys. Lett. A, 314:3, 222-231, 2003 

[12] A. Chen and J. Cao, "Almost periodic solution of shunting inhibitory CNNs with 
delays", Phys. Lett. A, 298:2-3, 161-170, 2002 



15 



[13] A. Chen and J. Cao, "Existence and attractivity of almost periodic solution for 
cellular neural networks with distributed delays and variable coefficients" , Appl. 
Math. Comp., 134:1, 125-140, 2003 

[14] C. He, "Almost periodic differential equations", Higher Education Publishing 
House, Beijing, 1992 



16 



